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ABSTRACT 

We analyze the color gradients (CGs) of ~ 50 000 nearby Sloan Digital Sky Sur- 
vey (SDSS) galaxies estimated by their photometrical parameters (Sersic index, total 
magnitude, and effective radius). From synthetic spectral models based on a simpli- 
fied star formation recipe, we derive the mean spectral properties, and explain the 
observed radial trends of the color as gradients of the stellar population age and 
metallicity. Color gradients have been correlated with color, luminosity, size, velocity 
dispersion, and stellar mass. Distinct behaviours are found for early- and late-type 
galaxies (ETGs and LTGs), pointing to slightly different physical processes at work 
qq \ in different morphological types and at different mass scales. 

In particular, the most massive ETGs (M» >10 n M©) have shallow (even flat) 
CGs in correspondence of shallow (negative) metallicity gradients. In the stellar mass 
range io 10 - 3 ~ 10 - 5 <^ M* <C, 10 11 M©, the metallicity gradients reach their minimum of 
f^*) | ~ —0.5 dex -1 . At M* ~ io 10 - 3 ~ 10 - 5 M©, color and metallicity gradient slopes suddenly 

change. They turn out to anti-correlate with the mass, becoming highly positive at 
the very low masses, the transition from negative to positive occurring at M* ~ 
10 9-9 ' 5 M©. These correlations are mirrored by similar trends of CGs with the effective 
radius and the velocity dispersion. We have also found that age gradients anti-correlate 
I— i ' with metallicity gradients, as predicted by hierarchical cosmological simulations for 

ETGs. On the other side, LTGs have gradients which systematically decrease with 
mass (and are always more negative than in ETGs), consistently with the expectation 
from gas infall and supernovae feedback scenarios. 

Metallicity is found to be the main driver of the trend of color gradients, especially 
for LTGs, but age gradients are not negligible and seem to play a significant role too. 
Owing to the large dataset, we have been able to highlight that older galaxies have 
systematically shallower age and metallicity gradients than younger ones. 

The emerging picture is qualitatively consistent with the predictions from hydro- 
dynamical and chemo-dynamical simulations. In particular, our results for high-mass 
galaxies are in perfect agreement with predictions based on the merging scenario, while 
the evolution of LTGs and younger and less massive ETGs seems to be mainly driven 
by infall and SN feedback. 

Key words: dark matter - galaxies : evolution - galaxies : galaxies : general - 
galaxies : elliptical and lenticular, cD. 
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1 INTRODUCTION 
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Galaxy formation is a complicated matter which has not 
yet come to a complete and coherent understanding. The 
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standard cosmological paradigm, the so-called ACDM, pre- 
dicts that dark matter (DM) haloes evolve hierarchically 
since early epochs, with smaller uni t s merging into massive 
structures ifeauffmann et al. 19931 . de Lucia et al. 200a . 
iRuszkowski fc Springelll2009h . which find positive evidences 
in the strong size evolut i on in massive galaxies since z ~ 3 
iGlazebrook et all |2004 iDaddi et al] 120051 . iTruiillo et aTl 
2006). This scenario is at variance with the evidences that 
today's high mass galaxies formed most of their stars in 
earlier epochs and over a shorter time interval than low- 
mass ones. This "downsizing" scheme (Cowic et al l 1 19961 ; 
Gallazzi et al. 1120051 ; [jimenez et al.ll2005): iNelan et al.l | 2005 : 
Thomas et alj 120051 : iTreu et al. 1 120051 ; iBundv et al] 2006 ; 



Cimatti et alj|200d; IPannella et al.ll2006l ; IPanter et al.ll2007 
Fontanot et al.ll200d ) seems at odd with the ACDM hierar- 



chical growth and seems rather to support a "monolithic- 
like" formation scheme (see e.g. Chiosi & Carraro 2002) 
which broadly predicts an inside-out formation of stars after 
an early dissipative collapse. 

However, it is increasingly clearer that the downsizing 
is compatible with the hierarchical model if the feedback 
processes in the mass accretion history of the ACDM 
dark h aloes are taken i nto a ccount (see e.g., iNeistein et al] 
20061. ICattaneo et alj l200St Conroy & Wechsler 2008, 



Cattaneo et alj|20ld). On the one hand, the gas cooling and 



the s hock heating ( Dekel fc Birnboiml l2006. Cat taneo et al] 
2008) are the main drivers of the star formation (SF) 
activity during the hierarchical growth of the DM haloes; 
on the other side supernovae (SNe), active galactic nuclei 
(AGN), merging, harassment, stran gulation etc. , gen- 



erally inhibit the stellar formation (iDekel fc Silk 19861: 



Recchi et al. 200ll; Pipino et alj 120021; Scannapieco et alj 
20061: iDi Matteo et all 120051: Ide Lucia et al. I 120061: 



Cattaneo et al. 2006: iKavirai et al.l 2007; ISchawinski et al 
20071; ICattaneo et al]|200Sl; lAntonuccio-Delogu fc Silkll2008l 
Khalatyan et al] 120081; [Romeo et al] |200S| ; iTortora et al 
2009bl ; IWeinmann et al. 1|2009T ). 



The different processes might rule the star formation 
at the global galaxy scale, or act at sub-galactic scales (e.g. 
the nuclear regions vs outskirts) such that they are expected 
to introduce a gradient of the main stellar properties with 
the radius that shall leave observational signatures in galaxy 
colours. 

This paper is motivated by the fact that color gradi- 
ents (CGs) are efficient markers of the stellar properties 
variations within galaxies, in particular as they mirror the 
gradients of star ages and metallicities (|Sandagelll972l ). al- 
though it is not yet fully clear whether met allicity is the 
main driver of the CGs in normal spheroidal dSaglia et al. I 
200d; ITamura fc Ohtal l2000l: lLa Barbera et alii" 



Kobavashil 12004 ISpolaor et al] 120091: 



and l ate-type systems (IMacArthur et al 



2002 



Rawle et alj 120091) 
2004 iTavlor et al 



2005 ). or age plays also an sig ni ficant role dSaglia et al. 
200ol; lLaBarbera et all 120031; iMac Arthur et alj I200I 



Spolaor et al]|200d ; iRawle et alj|2009l ). 

CGs are primarily a tool to discriminate the two broad 
formation scenarios (monolithic vs hierarchical), but more 
importantly they provide a deeper insight on the different 
mechanisms ruling the galaxy evolution. As these mecha- 
nisms depend on the galaxy mass scale, the widest mass 
(and luminosity) observational baseline is needed to remark 
the relative effectiveness of the different physical processes 



and their correlation with the observed population gradi- 
ents. 

The observational picture is so complex to justify a 
schematic review of the evidences accumulated so far. The 
facts are as follows. 

1) On average, nearb y elliptical and spiral galax- 
ies ar e bluer outwards ( Franx. Illingworth fc Heckrnanl 
I l98d. iPeletier et alj Il990al lbl iBalcells fc Peletierl 



outskirt s (IVader et 
IChabov" 



1994, 



ITamura fc OhtaT 20031). wh i le dwarfs show mainly redder 
aJjll983. iKormendv fc Diorgovskilll989l. 



1994 iTullv et al.Hl996l , iJansen et al]|200Ch . Re 



cently, deeper investigations of later-type galaxies have 
been done and have highlighted the presence of non- 
monotoni c colour and stellar population gradients in these 
syste ms 



ic colour and stellar population gradients m tnese 
jMacArthur et al] 12004 iBakos. Truiillo fc Pohlenl 



I2008L IMacArthur et al]|2009l. iMartmez-Serrano et al] 120091 
Sanchcz-Blaz quez et al. 1120091 ) . 

2) Spectral line indices (mainly measured in early-type 
galaxies, ETGs) are somehow a more efficient tracer of the 
stellar population properties and generally found to change 
with the radius like CGs |Kobavashi fc Arimotol Il99l 
Kunt schner et alj|2006l . IMacArthur et alj|200d . I Rawle et al 
2009h . 



3) Interpreting CGs in terms of metallicity gradients, 
in high-mass galaxies typically dlogZ/dlogR ~ —0.3, 
which is shallower than predicted by simulations of dissipa- 
tive "monolithic" coll apse (where —0.5 or a steeper value 
is expected, see e.g. Larson _ 1974 19751 . ICarlberj Il984 
lArimoto fc Yoshiilll987l . iKobavashilbOO^ ee also below). 

4) There are contradictory evidences of a correlations 
of CGs with galaxy mass and luminosity. Earlier studies 
have claimed a weak correlation with the p hysical proper- 
ties of galaxies (e.g., mass, l uminosity, etc., Peletier et alj 



1990al. 



Davies et alj iKobavashi fc Arimotol Il999l . 

Tamura fc Ohtal 120031 ), at variance with the typical mono- 
lithic collapse predictions. Recently, a stro nger correla- 
tion with mass has started to emerge (e.g. iForbes et al] 
2005), pointing to a metalli city gradient decreasing with the 
mass for low-mass galaxies (|Spolaor et al ] |200d . lRawle~ et alj 
2009), accordingly with the monolithic scenario. Instead, 
high-mass galaxies show too shallow gradients to match the 
predictions of the monolithic collapse, but compatible with 
galaxy merging. 

This wealth of observational evidences must be con- 
fronted with the model predictions coming from the different 
galaxy formation scenarios. 

1) Steep gradients are expected when stars form dur- 
ing strong dissipative (monolithic) collapses in deep poten- 
tial wells of galaxy cores where the gas is more efficiently 
retained, with a consequent longer star formation activ- 
ity and a longer chemical enrichment in the inner than 
in outer regions (negative metallicity gradients). On top of 
that, the delayed onset of winds from supernovae, causing 
a further metal supply in the central regions, would con- 
trib ute to reinforce th e steepness of these gradients (see, 
e.g. iPipino et al]|2008l ). As these processes are regulated by 
the galaxy potential depth, which is somehow related to the 
galaxy mass (and luminosity) Q, high-mass galaxies are ex- 



1 It remains to see whether there might be a correlation with 
the the DM fraction of the systems which is also a function 
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pected to have metallicity gradients which are steeper than 
lower-mass on e s. The latter seem to have almost no gradients 
dGibsonl 1 1991 IChiosi fc Carraro I l200l IKawata fc GibsonI 
2003), although the results here are based on limited galaxy 
samples. 

2) Within the hierarchical picture, merging induce a 
meshing of stellar populations, and lead to a more uniform 
metallicity distr ibution and to shallower CGs l|Whitelll98d 
IKobav ashi 2004) . Strong jets from powerful AGNs can also 
quench the star formation on the global gal axy scale and 
flatten colors gradients in the host systems |Tortora et al.l 
l2009bh . 

3) Metallicity and color gradients seem also to de- 
pend o n the efficiency of th e dissipative processes in dark 
haloes (|Hopkins et al.ll2009a), with just a weak dependence 
on the remnant mass llBekki fc Shioval Il999l . IKobavashil 
l2004l . IHopkins et ai1l2008bh . in a way more similar of the 
observed gradients. However, while merging is crucial to 
smear out the co lor radial variation in high-mass galaxies 
llKobayashill2004h, it seems unimportant for low- mass sys- 
tems ( de Lucia et al. II2006I . ICattaneo et a.1.1 l20Qgf ) . Instead, 
the energy from stellar winds and supernovae and the ef- 
fect of dis sipative collapse might induce even positive steep 
gradients (|Mori et al.lll997T l. 



CGs are the most direct observables to investigate the 
effect of the physical processes (such as merging, AGN, SNe, 
stellar feedback) , which drive the galaxy evolution as a func- 
tion of the main galaxy parameters: luminosity, mass and 
central velocity dispersion. In this paper we analyze opti- 
cal CGs in about 50 000 SDSS local galaxies, spanning a 
wide range of luminosities and masses. CGs are obtained 
in an indirect way from the photometrical parameters of 
the individual galaxies (effective radius, Sersic index, total 
magnitude). Using synthetic spectral models, we compute 
age and metallicity gradients for early- and late-type sys- 
tems, and analyze the trends with mass. The sample allows 
us to investigate the distribution of the CGs over an un- 
precedented baseline of galaxy sizes, luminosities, velocity 
dispersions, and stellar masses. The observed trends are in- 
terpreted by means of quite different physical phenomena 
at the various mass scales. This approach, while exposed to 
the uncertainties on the structural parameters (e.g. Sersic 
index, n, and effective radius, -R e ff) has the advantage of 
dealing with large statistics; furthermore it is exportable to 
higher redshift galaxy surveys which are generally limited 
to rest-frame visual bands. In this respect we pay particular 
care in the check of the consistency of our results with inde- 
pendent analyses. We will mainly concentrate on ETGs, for 
which a wide collection of results (having an homogenized 
gradient definition) is available from literature. On the other 
hand, due to the uncertainties of the structural parameters 
of late- type galaxies (LTGs), we will take the results on this 
galaxy sample with the right caution, being well aware that 
our findings might be only a benchmark for more accurate 
analyses. 



In Sj2]we present the data sample and spectral models. 
In £[3] we show the main results of our analysis, discussing 
CGs as a function of structural parameters, stellar masses, 
and stellar properties derived from the fitting procedure. In 
the same Section, age and metallicity gradients are shown, 
and discussed within the galaxy formation scenarios in i}?] 
We finally draw some conclusions in ij5] 

In the following, we use a cosmological model 

with (fi m , Q A , h) = (0.3, 0.7, .7), w here h = 

iio/100 km s^Mpc -1 (|Spergel et all 120031 . 12001 . cor- 
responding to a Universe age of t un iv = 13.5 Gyr. 



2 DATA 

2.1 Galaxy sample 

Our database consists of 50 000 low redshift (0.0033 ^ 
2 < 0.05) galaxies in the NYU Value-Added Galaxy Cat- 
alog extracted from SDSS DR4 jBlanton et al.ll2005bl . here- 
after B05fl B05 have performed a set of quality checks 
contemplating the analysis of large and complex galaxies, 
incorrectly deblended in the SDSS pipeline, and of the 
star/galaxy separation, with a number of eyeball inspec- 
tions on objects in the catalogue. The catalog includes 
Petrosian magnitudes and structural parameters in ugriz 
bands such as the concentrations C = R90/R50, where Rgo 
and R50 are the Petrosian radii enclosing 90% and 50% of 
total luminosity. Effective radius, R c s, and the Sersic in- 
dices n are taken from B05, where, a ID seeing-convolved 
Sersic profile is fitted to galaxies. The adopted redshift 
range sets the completeness level at a Petrosian magni- 
tude r ~ 18 and at a surface brightness within R50 pso — 
24.5 mag arcsec -2 . This luminosity limited sample still suf- 
fers from an incompleteness at low surface brightness levels, 
e.g. in the regime of dwarf galaxies with M, <i, iq 8 - 3 ~ 8 - 5 JVf ( 



ilBlanton et"afl l2005al . iBaldrv. Glazebrook fc Driver! 120081 . 
iLi fc White! 20091 ). Moreover, the incompleteness, almost 
negligible for more massive galaxies, comes back again for 
the most massive ones (M* > 10 Mq) as the sample is bi- 
ased against large objects at lower redshifts i n consequence 
of a b ias in the SDSS photometric pipeline (|Blanton et al.l 
l2005al ). We address the systematics of the SDSS sample 
more in detail in App. |A"T1 

Petrosian magni tudes are reduced to 2 = (K-term; 
iBlanton et alj 2003a!), corrected for galacti c extinction us- 
ing the Schlegel. Finkbeiner fc Davisl dl998f) dust map , con- 
verted to our adopted cosmology ( Blanton et al.ll2003bh , and 
further corrected for the missing flux of Petrosian magni- 
tudes. The latter correction is a lmost unchanged using the 
relation in iGraham et al] (|2005l ). 

The sample spans a wide range of luminosities, masses, 
and colors, including galaxies lying both in the red-sequence 
(RS) and in the blue cloud (BC). In particular, we have 
sorted out ETGs by adopting the following two criteria: 

• the Sersic index satisfies the condition 2.5 ^ n ^ 5.5; 



• the concentration index C > 2.6 ([Sh imasaku "et al.l 
l200ll . Padmanabha n et al. l l2004ICardone et al . 2008)0. 



of the galaxy stellar mass (see e.g., |Napolitano et al. I |2005l. 



Cappellari et al. 1 12003 , IConrov fc Wechsle j |200S|. iTortora et alj 
2009al) . 



The 



catalog 



available 



at: 



http : //sdss .physics .nyu. edu/vagc/lowz .html 

a In order to retain the low mass red systems, no lower cut has 
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Figure 1. Gradients of the g — i colour, as a function of various observed and derived quantities. Medians of binned CGs are plotted 
together with the 25 — 75% quantiles, shown as error bars. Top-left. CG as a function of color, black open and gray squares are for colors 
at i?i and i?2, respectively. Top-right. CG as a function of r-band magnitude. Middle-left. CG as a function of logarithm of ctq, which is 
defined as the veloci ty dispersion within a c i rcular aperture of radius 0.1-R e ff (using r-band i? c fj) from the SDSS velocity dispersion <r ap , 
using the relation in I j0rgensen et "aL I jl995l , ll996l) . Middle-right. Gradient as a function of total stellar mass M t (assuming a Chabricr 
IMF), which is the output of our stellar population analysis fitting the total colors. Bottom-left. CG as a function of logarithmic of 
r-band R c ff (the results does not depend on the band). Bottom-right. CG as a function of i-band Scrsic index n; the lines are the the 
median trends for the Sersic index in the other bands. 
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The final ETG sample consists of 10 508 galaxies. Of the re- 
maining entries of our database, 27 813 are LTGs, defined as 
objects with C 2.6, n ^ 2.5, and oo ^ 150 km/s, and 9 359 
are intermediate morphology galaxies, with C ^ 2.6 (> 2.6) 
and n > 2.5 2.5), hereafter named Intermediate- type 
galaxies (ITGsjS Further details on the adopted morphol- 
ogy criteria are given in the App. IA2I In this paper we will 
consider ETGs and LTGs only, with a main focus on ETGs. 



2.2 Spectral models and gradients 

The stellar population analysis of the selected sample is 
based on a set of "single burst" synthetic st e llar s pectra 
using the prescription of iBruzual fc Chariot I |2003l . here- 
after BC03). These models are widely used in literature 
both for ellipt i cals a n d spirals (e.g . . iTrager et ail 1200(1 
Bell fc de Jongl l200ll . iGanda et ail 120071 iPeletier et al l 
20071 ). but we are aware that, in some cases, e.g. in 



later-type systems wi th very protracted star formations 
l|Kennicutt et al.lll983l). they may fail to repro duce correctly 
the galaxy properties ( MacArthur et al1l2009l ). A finer anal- 
ysis would require a more realistic SF history for LTGs (e.g., 
an exponentially decaying or a constant SF), which is be- 
yond the scope o f this pape r. 

We adopt a IChabrier I (|200ll , l2002i , l2003h initial mass 
function (IMF), noting that the choice of a lSalpeter I (|l955l ) 
IMF would affect the total stellar mass only, shifting log M* 
upward by ~ 0.25 and leaving a ll the other stellar p opula- 
tion parameters unchanged (e.g. iTortora et all 12009a ), Our 
choice of the BC03 models is driven by their high versatility 
and ability to span the stellar parameter spa ce, e.g. metal- 
licitie s and ages. There are other options (|Rettura et al, 
I2006L iKannappan fc Gawiserl 120071 . ITortora et al.1 12009a 
for a comparison of the results from di fferent prescrip- 
tions), but none is found to be bias free (iMarastonl 20051. 
van der Wei et aiT 2006, Maraston et al. 2009, IConrov et al.l 



20091 . l2010al lbir 

In our model procedures we adopt a metallicity ranging 
within 0.0001 — 0.05, and the Universe age, t un iv (in the as- 
sumed cosmology), as an upper limit for the galaxy ages. 
Synthetic colors have been obtained after convolving the 
spectra with the corresponding SDSS filter bandpass func- 
tions. To improve the sensitivity of the results to the small 
differences in both Z and ages, we have adopted a "mesh re- 
finement" procedure which interpolates the synthetic mod- 
els over the initial grid of equally spaced bins in the loga- 
rithm of Z and age. To estimate the galaxy age, Z, and the 
stell ar mass-to-light ratio , T„, we use a % 2 minimization fit 
(see lTortora et a l. 2009a b for details) to the observed colors 
(u — g, g — r, g — i, and g — z). It is known that the use of just 
the optical broad bands is severely prone to the well known 
age-metallicity degeneracy. Ho wever, following previous in- 
dications (e.g. IWu et al.l I2005T ) we have verified in App. B 
that the use of near infrared (NIR) photometry reduces the 



been applied to the velocity dispersion. However, very low val- 
ues (ao < 70 km/s) may be prone to large systematics, due to 
the signal-to-noise ratio and instrumental resolution of the SDSS 
spectra. 

4 Note that 2120 galaxies are cut away by the conservative upper 
limit on n, dictated by the fitting Sersic code limit of n = 6. 



overall uncertainties but does not modify the estimates of 
the intrinsic population parameters based on optical bands 
only. We shall return on this issue later on. 

We have used the structural parameters given by B05 
to derive the color profile (X — Y)(R) of each galaxy as 
the differences between the (logarithmic) surface brightness 
measurements in the two bands, X and Y. The CG is de- 
fined as the angular coefficient of the relation X — Y vs 
6{X - Y) 



log R/Rea, V; 



measured in mag/dex 



6\og(R/R aB ) 

(omitted in the following unless needed for clarity). The fit 
of each color profile is performed in the range Ri = 0.1-R e H ^ 
R ^ R2 = R e fi- We have verified that a choice of different 
radial ranges for the fit leaves the global trends unaffected. 
In particular, by setting R2 = 2i? c ff the gradients changed 
by <C, 0.01 — 0.02 mag, where R e s is the r-band eff ective ra- 
dius {Peletier et al.lll990bl . ILa Barbera et al.ll2005l ). Slightly 
di fferent definitio ns are also used elsewhere (see discussion 
in lLiu et afll2009h . 

By definition, a positive CG, \?x-y > 0, means that 
a galaxy is redder as R increases, while it is bluer outward 
for a negative gradient. The fit of synthetic colors is per- 
formed on the colors at Ri and R2 and on the total inte- 
grated colors. Following lLa Barbera et all {2005) we define 
the stellar parameters gradients as V age = log[age2/agei] 
and Vz = log[Z2/Zi], where (agd, Zi) with i = 1, 2 are the 
estimated age and metallicity at Ri and R2, respectively. 

A full test of the reliability of our modelling technique, 
the presence of spuriously generated correlations, and the 
contribution from internal dust extinction are addressed in 

App.m 



3 RESULTS 

We show how color and stellar population gradients depend 
on structural parameters and mass, and compare the re- 
sults with the predictions from galaxy formation scenarios. 
The median trends, e.g. in bins of mass, velocity dispersion, 
colors etc., will be represented along with the 25-75th per- 
centiles in each bin. 

The average behaviour of galaxies is analyzed in the 
parameter space, and the galaxy properties discussed after 
splitting the sample in the two morphological classes, ETGs 
and LTGs. We concentrate mainly on the analysis of ETGs 
and compare our findings with a wide collection of indepen- 
dent data and simulations. 



3.1 CGs as a function of structural parameters, 
luminosity and mass 

Fig.[T]presents the gradients V g _i as a function of integrated 
colors, structural parameters, r-band magnitude, central ve- 
locity dispersion, and stellar mass. We adopted the g — i 
color as a reference having varified that the use of differ- 
ent colors does not affect our conclusions. As shown in the 
top-left panel of the figure, CGs are, on average, positive 
for the galaxies with bluer colors, g — i < 0.5, while they 
reach increasingly negative values f or redder galaxi es (e.g. 
iPark fc Choill2005l ; iLee et al.ll2007l : fSuh et aJj|2uTch . Since 
g — i~ 0.5 marks the transition from the RS to the BC, this 
result allows us to broadly associate negative and positive 
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Figure 2. V g _i as a function of effective radius (left panel), velocity dispersion (middle panel) and stellar mass (right panel) for ETGs 
(red symbols) and LTGs (blue symbols). 



Table 1. Slopes of the correlation between CGs (V u — g , V g — r , 
V g _i, Vg_ z ), Vage and Vz vs log Co and logAf* for ETGs and 
LTGs. The errors on slopes are computed as la uncertainty by 
a bootstrap method. * For the LTG sample, splitting galaxies by 
masses smaller /larger than logAf* ~ 10.3 would result in slopes 
of -0.07 ± 0.01 and 0.26 ± 0.25, respectively. 



Correlation 







log cr < 2.2 


log cr > 2.2 






- log a 


-0.60 ± 0.03 


0.40 ± 0.05 


0.02 ± 0.05 


V E -r 


- log a a 


-0.27 ± 0.01 


0.19 ± 0.01 


-0.02 ± 0.02 


Vg-i 


- log CTq 


-0.34 ± 0.02 


0.25 ± 0.04 


-0.03 ± 0.02 




log CTQ 


-0.40 ± 0.02 


0.17 ± 0.04 


-0.09 ± 0.03 


Vago 


- log <T 


-0.10 ± 0.03 


0.04 ± 0.17 


0.14 ± 0.02 


V Z - 


log CTQ 


-0.86 ± 0.07 


1.02 ± 0.24 


0.06 ± 0.10 



Vg 
Vg 
Vg 

V a 



, - log JVf, 
* - log M, 

- log M, 
, - log M, 

- log M„ * 
■ log M, 



logM, < 10.3 
-0.43 ± 0.04 
-0.24 ± 0.02 
-0.34 ± 0.02 
-0.36 ± 0.03 
-0.18 ± 0.06 
-0.46 ± 0.06 



logM, > 10.3 
0.12 ± 0.02 
0.05 ± 0.01 
0.07 ± 0.01 
0.02 ± 0.01 
0.24 ± 0.03 
0.10 ± 0.04 



-0.30 ± 0.01 
-0.13 ± 0.01 
-0.18 ± 0.01 
-0.18 ± 0.01 
-0.01 ± 0.01 
-0.45 ± 0.01 



CGs to the RS and BC, respectively. Furthermore, galaxies 
with intermediate colors (g — i ~ 0.5 — 0.7 or equivalently 
g — r ~ 0.2 — 0.6) have almost null (or slightly negative) 
gradients. The same trend holds by considering central and 
outer colors (at R = Ri and i?2). We will return to this issue 
later on. 

Tight trends are found when plotting the gradients 
against the r-band total magnitude and stellar mass, while 
the dependence on ctq is weaker. In the latter case, gradients 
look negative almost everywhere with the less massive galax- 
ies (log (Jo < 2 km/s) having the smallest V g _i (~ —0.2). 
Although our sample does not include very massive galaxies 
(M, ^lO 11 ' 5 Mq), our results seem to show a stable trend at 
the large mass scales, pointing to even shallower gradients. 
We note here that the smaller range shown by V g -i as a 
function of tro in the central panel (left column) of Fig. [I] is 
mainly due to the mix of the ETGs and LTGs. They have 
different trends with small cro in the low mass regime (see 
also i]3.2[) . thus producing a dilution of the average gradients 
with respect, e.g., the larger median V g _i (~ 0-^0.2) found 
for the low stellar mass (right column). 

The correlations with luminosity and mass show that, 
starting from the bright/massive end, the CGs become 
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Figure 3. CG V g — r as a function of r-band m agnitude for ETG 
sample. Superimposed are the results from lLa Barbera et all 
(2005) for high-richness (open boxes), low-richness (filled boxes), 
and all (crosses) the clusters. The blue and red a rrows are the me - 
dian V g — r for the massive g alaxies analyzed in lWu et al. I fe005l l 



and lLa Barbera et al.l 1 120091 ) respectively; see the text for details. 



steeper, with the steepest negative gradients (~ —0.2) cor- 
responding to r ~ —20 mag and M» - 10 la3 M Q . From 
there, bluer and later-type systems begin to dominate, and 
the gradients invert the trend with luminosity and mass, 
becoming positive at r ~ —17.5 mag and M* ~ 1O 8 ' 7 M0. 
The lowest mass systems are those with the highest, positive 
gradients (~ 0.2). The mass scale at which the CG trend in- 
verts is in remarkable agreement with the typical mass scale 
break for "bright" and "ordinary" galaxies (Capaccioli et al. 
1992, see also jHT^l) or for star-forming and pass ive systems 
i|Kauffmann et alll2003l . iDrorv fc Alvarez! 120081 ) . 

The above trend holds also for the correlation with 
R a H and, less strongly, with the Sersic index, n, as 
expected in view of the tight correlation with stel- 
lar mass and luminosity llCaon. Capaccioli fc D'Onofriol 
| l99St IPrugniel fc Sirnienl I1997I. iGraham fc Guzman! 120031 . 
Ishen et al.li2003l . iTruiillo et al. |[2004 iTortora et al.ll2009aF ). 
Note that the minimum at n ~ 2 for the trend of the gradi- 
ent (Fig. [1] bottom right) roughly coincides with the value 
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Figure 4. Age and metallicity gradients as a function of velocity dispersion (top panels) and stellar mass (bottom panels). Solid and 
dashed lines are for galaxies with central stellar populations older and younger than 6 Gyr, respectively. Symbols are as in Fig. [2] 



used here to set the separation between ETG and LTG sys- 
tems, suggesting differences in the CG properties. 

3.2 Gradients in morphologically selected galaxies 

Following the dependence of V g _i vs n just discussed, ETG 
and LTG samples will be now considered separately. In Fig. [2] 
we plot V g _i versus R e g, cro, and M t for ETGs (red symbols) 
and LTGs (blue symbols). The slopes of the linear best-fit 
relations of the CGs versus cro and M* are shown in Table 

in 

First, ETGs are on average br i ghter and more mas- 
sive than LTGs (|Kauffmann et alTI I2003T ) and have less 
dispersed CGs. As a general result, large R e g galaxies 
are found to have lower gradient s than more compact 
systems (e.g. I Hopkins et al.l l2009al ). However, for LTGs, 
Vg_i monotonically decreases with R B g and becomes flat 
at log R c ff J>0.5, while for ETG it slightly increases for 
logi? ff ;>0.3, and decreases for smaller R c h. Similar trends 
are observed for the gradients as a function of the stellar 
mass. LTGs show a monotonic decreas ing t rend (in qualita - 
tive agreement with iKim fc Annlll99d and lLiu et al T l2009T l. 
while a U-shaped function is found for ETGs with the 
gradients definitely decreasing with the mass for M* <C, 



10 10 ' 3-10 ' 5 Mq, and a mildly increasing for larger mass val- 
ues. This mass scale roughly corresponds to a total luminos- 
ity of r ~ —20.5 mag, which is compatible with the typi- 
cal luminosity (and mass) scale for ETG dichotomy in the 
galaxy structural properties (e.g.. ICapaccioli et al . 1992a, 
ICapaccioli. Caon fc D'Onofriolll992bl. iGraham et alj|2003h . 
As discussed in lCapaccioli et al.l l|1992al ) this dichotomy may 
be related to the formation mechanism processes, and in 
particular to the dominance of the merging events for the 
"bright" sample, which is compatible with the shallower gra- 
dients that we observe for these systems (see, e.g., Fig. [21 
right panel). Nevertheless, for a fixed stellar mass we observe 
that ETGs gradients are, on average, larger than LTG ones. 

The same two - fold trend is shown for ETGs gradients 
as function of the velocity dispersion, while no trend with 
the mass is observed for LTGs gradients. 

As a consistency check, in Fig. [3] we superpose to our 
g — r grad ients as a function of the r— band luminosity the 
results of lLa Barbera et all |2005) obtained by analyzing 
a sample of low redshift luminous (r ^ —20 mag) cluster 
galaxies. Their value (V g _ r ~ —0.075) is in agreement with 
our findings over the range where the two studies overlap, 
but we were also able to identify a change of V g _ r with 
our larger luminosity baseline. The results from the massive 
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galaxies in IWu et alj (|2005l ) and lLa Barbera et all (|2009l ) 
are also shown to have V g _ r consistent with our results. 

Unfortunately, we have no information on the envi- 
ronment of the individual systems in our sample, to check 
the dependence of our results with the local density. There 
are many lines of evidences that ETGs in clusters (and in 
very dense environments in general) have sh allower gradients 
than systems in less den s e environments (Tamura fc Qhtal 
120001 . iTamura et al]|2000l . lLa Barbera et alj|2005l - see also 
Fig. O- However, according to the morphological segrega- 
tion, LTGs, which are expected to reside in the field or in 
the external regions of clusters, have indeed lower gradients 
than ETGs which are found to live mainly in higher density 
environments. 



3.3 Interpretation in terms of Z and age gradients 

Here we test whether the CGs can be explained in terms of 
stellar properties such as age and metallicity. As seen in 
we have obtained the gradients of metallicity and age distri- 
butions using two characteristic scales (0.17? e ff and li? c ff). 
The V agc and Vz are shown in Fig. U as a function of the 
velocity dispersion and the stellar mass, and the fitted slopes 
of such relations are quoted in Table [1] 

Some features survive the large scatter. For LTGs the 
Vag C is about zero both with 00 and mass, strongly sug- 
gesting that CGs should not depend on age gradients of the 
galaxy stellar population^. Actually, Vz is instead strongly 
dependent on A/*, with the lowest metallicity gradients 
(~ —1) at the largest masses, while the dependence on ao is 
very weak as the Vz ~ (—0.5, —0.6) for all 00 • 

ETG age gradients seem basically featureless when plot- 
ted against <tq where we find V agc > with typical me- 



dian values of V a 



0.2. On the other hand, the metal- 



licity gradients show some features with Vz decreasing for 
log ao < 2.2km/s, with a slope = —0.86, and increasing for 
larger velocity dispersion, with a slope of 1.02, reaching the 
shallowest values (~ —0.2) at log 00 J>2.4km/s and at the 
very low <tq (i.e. log (To < 1.8km/s) where Vz~ 0. This pecu- 
liar trend is mirrored by a similar dependence on the stellar 
masfl This two-fold behaviour is also significant when we 
plot Vagc as a function of the stellar mass: V age decreases at 
log M* < 10.3 — 10.5, and increases in more massive systems. 
Bright ETGs with r <; -20 mag (logM* £10) have a 



median V aEC = 0.144±° °°* (with scatter and Vz = 

-0 004 y™^ 11 =^»"ci _o 27)' where the error bars rep- 



' a g e — w - XT1 -0.003 

-0.472±°'°°i (with scatter + ' ' 1 



resent the la uncertainty and the scatters the 25-75 per- 
centiles of the sample distribution. However, we have found 



5 The scatter of the gradients around zero has been suggested to 
be a manifestation of possible differences in th e formation pro- 
cesses in these systems MacArthur et al. (2009). 

6 For ETGs, stellar mass and ao are tightly correlated, thus the 
correlation and the interpretation discussed here can be made 
when dissecting the similar trend as a function of velocity disper- 
sion. On the contrary, many low mass LTGs have very large ve- 
locity dispersion which contribute to have almost constant metal- 
licity gradients for logcro < 2.2km/s and increasing at larger ctq. 
But in this case the scatter is very large suggesting that the cen- 
tral velocity dispersion is not a representative parameter of such 
rotational velocity supported systems. 
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Figure 5. Metallicity vs age gradients for ETGs and LTGs. The 
symbols are the same as in Fig.[2]and|4] Metallicity vs age gradi- 
ents for ETGs are compared w ith predictions from gas-rich merg- 
ings in Hopkins et al. ( 2009a) (yellow shaded region) and cluster 
ETGs in lRawle etafl (20051 ) Cblack points). 



that V a gc and Vz strongly depend on galaxy age: in par- 
ticular, if we separate the sample in systems with older and 
younger than 6 Gyr central stars, we obtain different age and 
metallicity gradient trends as shown in Fig. [4] Older systems 
have gradients which are shallower than younger (in partic- 
ular for ETGs) and bracket the average trend of the whole 
samples. A partial contribution to the scatter is also given 
by the central metallicity: the lower-Z systems show gener- 
ally shallower gradients and viceversa. However, the effect 
is significantly smaller with respect to the central age. 

The median values of the older sample is V agc = 01q 004 



0.337ln nns (with scat- 



(with scatter Z ; 14 ) and Vz — — u.u«ji_ 005 
tcr 

-oil) while for the 

younger one we obtain V 



0.43; 



(with scatter toTo) and Vz = -0.67±%%& (with 
scatter +0 13)- Galaxies with old central ages are fully 
consistent with the resul ts of pure metallicity a nd mixed 
age+metallicity models in |La Barbera et alj (l2005h and with 
the best fitted values in La Barbera et al. I l|2009l ) (where 
galaxy ag es are in the range of 7.8 — 12.6 Gyr). The agree- 
ment with lRawle et al. I (|2009h is also good if we consider that 
their sample has a median central age of ~ 10 Gyr (calcu- 
lated within an 7? e ff/3 aperture), comparable to the old ob- 
jects of our sample: they find an age gradient of —0.02 ±0.06 
and metallicity gradients of —0.13 ± 0.04. Finally, from the 
anal ysis of 36 nearby ETGs with SDSS and 2MASS photom- 
etry, IWueTal] (|2005h have derived V agc = 0.02 ±0.04 (and 
scatter 0.25) and Vz = -0.25 ± 0.03 (and scatter 0.19), 
which turn to be qualitatively consistent again with our 
"old" sample. The inclusion of dust extinction leaves our 
general considerations unaffected, as we will show in App. 
IB3I The average trends of the metallicity and age gradients 
are almost unchanged, with variations going in the direction 
of a better match with other literature data. 

The massive ETGs, with young cores (< 6 Gyr) have 
very steep metallicity gradients, which are discrepant with 
respect to the literature results we have compared with. In 
App. [A"2l we have checked that this result does not depend 
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on the contamination of the ETG sample from late-type sys- 
tems. These latter, indeed, also show the same dependence 
of gradients on central age (see Fig. |4]). A clear example of 
systems with very steep gradients are the isolated galaxy 
NGC 821 which have a young central stellar population and 
a metallicity gradient of Vz = -0.72±0.04 or NGC 2865, for 
which is V z = -0.47±0.05 l|Proctor et al.ll2005l , lReda et all 
120071 ) . One reason why we have been able to pick up this dif- 
ference between centrally old and young systems is that most 
of the samples ana lyzed in literature are made up of clus - 
ter ellipticals (e.g. iMehlert et"afl 12001 iRawle et al.l |2009h 
with only few c ases of field galaxies ijOeando et al. 1 120051 . 
iReda et al.ll2007l ), while our sample is composed by galaxies 
living in various environments. On average, field galaxies are 
foun d to be younger and have a wider distribution of galaxy 
age jBernardi et al.lll998l iTrager et al]|2000l . iThomas et all 
120051 ). thus, a large fraction of our centrally young galaxies 
would be made of systems in the field or low density environ- 
ments with steeper age and metallicity gradients than those 
in cluster galaxies (due to the little interactions with the 
neighborhood, see Sect. IqT21 and lLa Barbera et all (|2005l )). 

Fig. [S] displays a relevant anti-correlation between age 
and metallicity gradients: for ETGs the correlation is clear 
with a slope of —0.94 ±0.02, while for LTGs it seems poorer 
but still significant, with a slope —0.48 ± 0.03 (see also Ta- 
ble [2]). This is an intrinsic property of the galaxy sample, 
which is consi stent with predictions from the galaxy merging 
scenario (e.g. iHopkins et al.ll2009a|) and with other observa- 
tional evidences (e.g. Rawle et al. 2009). It is not spuriously 
induced by the well known age-metallicity degeneracy (see 
App. B). In particular, galaxies with strong negative metal- 
licity gradients have positive age gradients (coherently with 
the merging simulation in the figure), while the occurrence 
of negative age gradients is associated to small mass systems 
with positive metallicity gradients. The negative correlation 
found for LTGs is also qu alitatively consistent wi th the lit- 
erature on spiral galaxies l|MacArthur et al]|2004l ). 

In Fig. [5] we finally show the correlations between gra- 
dients and stellar parameters at different galaxy radii, with 
their related slopes listed in Table [5] Metallicity and age 
gradients show a tight correlation with central metallicity 
and age, respectively. For the ETGs, the age gradients are 
generally different from zero, positive for small central ages, 
and null or negative for old systems. Metallicity gradients 
are negative for systems with older stellar populations at 
.Reff or with higher metallicity, while they are about zero 
for systems with old central ages or intermediate ages at 
ReB (4 — 7 Gyr). Noticeably, the very positive Vz and neg- 
ative Vage correspond to very young populations at R e g, 
i.e. recent star formation outside the galaxy cores. These 
latter systems seem c ompatible with d warf galaxies hosting 
expanding shells (e.g. iMori et ah 11.1991 see also ©. 

Similarly, LTGs have also positive age gradients in pres- 
ence of young cores, but in these cases we observe also large 
central metallicities (see Fig. [6] top row, third panel from 
left) and for these systems we also found the strongest metal- 
licity gradients (bottom raw, third panel of Fig. [6]). 

Overall, results in Figs. [4] and [6] are i n substantial 
agreement with the correlations reported in I Rawle et al.l 



|2009l ) and with the findings in hydrody namical simu- 
lations of both elliptical and disc g alaxies (|Hopkins et al.l 
l2009al . ISanchez-Blazquez et al.ll2009l ) . 



3.4 Comparison with literature data 

We now proceed to a more detailed comparison of our 
findings with a set of literature works which make use 
of a more sophisticated analysis, although usually asso- 
ciated to smaller samples. In particular, we concentrate 
on the ETG sample, taking the results for LTGs with 
the necessary caution, due to the larger degree of un- 
certainties arising from, e.g., 1) the simple assumption 
on the star formation recipe, 2) the assumption of no 
dust gradients (see also App. IB3|) . and 3) the prediction 
from theoretical simulations of the formation of such com- 
posite and complex systems which are still in their in- 
fancy jMacArthur et al]|2009l iMartmez-Serrano et al.ll2009l . 



ISanchez-Blazquez et al. 2009) 

Due to our minimal choice of the wavelength base- 
line and the possible biases on the stellar population es- 
timates caused by the use of the optical band only, the 
comparison of our results with literature data is critical 
to 1) ultimately check the robustness of our results on 
the regions of the parameter space overlapping with in- 
dependent analysis and 2) assess the gain in the physical 
information allowed by our analysis. In Fig. [7] we sum- 
marize the metallicity and age gradients for ETGs as a 
function of central velocity dispersion, and compare these 
with the observed gradients f rom a collection of literature 
studi es dMehlert et all 120031. iProctoij 120031. lOgando et al.l 



20051. IReda et all 



jgOOTj 



Sanchcz-Blazqucz et al. 



2007 



Koleva et all l2009al lbl. ISpolaor et al.l 120091 . iRawle et all 
20091 ). In the next section we will expand this comparison to 



the predictions from cosmological simulations. When consid- 
ering massive ETGs, our Vz are systematically steeper than 
literature data (see top- left of Fig. [7|), although showing on 
average a lower scatter than dwarf galaxies. 

As seen in Fig. [()] and also illustrated in Fig. [7] (top- 
right panel), a significant contribution to the scatter of the 
metallicity gradients is linked to the age of galaxy cores, 
with older galaxies having, on average, shallower gradi- 
ents than younger systems. When considering objects with 
agei > 6 Gyr, the agreement with the other studies (gener- 
ally dealing with old systems; see H3.3P is remarkably good. 
The same conclusion is reached when considering the V age 



7 In particular, if we limit our galaxies sample to luminous sys- 
tems (r < —20 mag) with Z\ ^ 0.008 and agei ^ 2 Gyr and per- 
form the fit with respect to the central quantities within R e ff/3 
we find that the slope of the Vz — log Z\ relation is —0.58 ± 0.04 
in good agreemen t with the values in the range [—0.6, —0.7] from 
I Rawle et alj (I2009T) . For V ag e — log age\ relation we find a slope of 
— 0.90 ± 0.06, which d isagree with the shallower values reported 
in IRawle et al. I J2009[) . A similar disagreement is found for the 
Vz —log agei and V age — log Z\ correlations, where we find slopes 
0.67 ± 0.05 and 0.50 ± 0.03 respectively, which are still not per- 
fectly consistent with their results. These might be related to the 
small sample adopted, to the selected environment (i.e. cluster 
cores), to differences in population markers (i.e. line strengths 
rather than colors) and the use of synthetic spectral models with 
variable o-enhancement over-abundance. 
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Table 2. Slopes of the linear correlations between stellar population gradients Vz, V agc , Z and age for ETGs and LTGs. IN and OUT 
indicate gradients correlated with inner (at Ft = Ri) and outer (at R = R2) metallicity and age. In the last row the result for correlation 
between Vz and V age are reported. The fitting procedure is the same as that in Table [T] W; here the linear correlation is a 0-th order 
approximation. 



Correlation 7ETG 7LTG 
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OUT 
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- log age 
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0.75 ±0.03 
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-logZM 
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V z - 


log age 
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-0.20 ±0.06 
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Figure 6. Correlations of gradients and fitted stellar parameters. The symbols are the same as in Figs. [2] [4] and [5] Top panels. From 
left, age gradient as a function of age\, age2, Z\ and Z2. Bottom panels. From left, metallicity gradient as a function of age\, age2, Z\ 
and Z2. 



(Fig. bottom-right), where the steeper (positive) gradi- 
ents are shown by the younger systems. In App. lB3l we show 
that shallower age gradients might be found if the dust ex- 
tinction is accounted in the stellar models for the massive 
ETGs, while metallicity gradients remain almost unchanged. 

The dependence of the gradients on the galaxy age is 
consis tent w ith results in t he dwa rf regime bv lSpolaor et alj 
(|2009h and iKoleva et all |2009al lbh. which find shallower 
(and tightly correlated with mass) and steeper gradients in 
old and young systems, respectively. 

As a final remark, galaxies with very low velocity dis- 
persion (log (Jo < 1.85 km/s, see the vertical dashed gray line 
in Figs. [7] and [8]) could be affected by some systematics due 
to the low signal-to-noise ratio and instrumental resolution 
of SDSS spectra. Notwithstanding this caution, our recov- 
ered gradients rep roduce fairly well th e literature results for 
dwarf ETGs (e.g. ISpolaor et all 120091 ). Also, the indication 
of shallower gradients at the lowest 00 seems robust, due to 
the clear turnover at log 00 ~ 2.2 km/s and the inversion in 
the range 1.85 <C, log 00 ^ 2.2 km/s. 



4 DISCUSSION 

Our results seem to support the idea that the metallicity 
trend versus the stellar mass for LTGs is mainly driven by 
the interplay of gas inflow and winds from supernovae and 
evolved stars <|Matteuccilll994l . lKobavashil2"ool . iPipino et al.l 
2008). These processes tend to increase the central metallic- 
ity and prevent the enrichment of the outer regions. There- 
fore more massive systems have on average larger central 
metallicities which correspond to steeper negative gradients. 

Low mass ETGs (M» < lO 10 ' 3-10 ' 5 M Q , see Tabled 
show a similar correlation of the metallicity gradient with 
stellar mass (Fig. 2]), which suggests that they might experi- 
ence the effect of infall/SN feedback as LTGs. For instance, 
low mass ETGs in Fig. [8] are consistent with SN feedback 
(both soft an d strong ) , as predicted in chemo-dynamical 
simulations in iKawatal l^OOlT ). On average, weaker SN feed- 
back gives gradients in agreement with ours at intermediate 
log (Jo (~ 2.2 km/s), while the stronger feedback recipe seems 
to reproduce better the low log 00 side of the correlation. A 
recipe including the SN feedback with a varying power as 
a function of mass l|Dekel fc Silklll986l . iDekel fc Birnboiml 
2006) would allow a fairly good match of the observed de- 
creasing trend of the gradients for the less massive galaxies 
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Figure 7. Metallicity gradients for ETGs as a function of central velocity dispersion compared with some literature data. Galaxies with 
velocity disper sion < 70 km/s, at the left of the vertical dashed gray line, are somehow uncertain. Top-left. Me tallicity gradients are 
compared with iRawle et aT I J2009l) (R09) as black poin ts (including the results for Shapley and A33 89 clusters) , [Mehlert et all i2003ri 
(M03) as blue triangles. ISanchez-BIazq uez et al. I ll2007f) (S B07') as purple circles. FSpolaor et al-l d 200Sl') ( S09) as green stars. lKolevaet**aD 
l2009allb | ) (K0 9) as light green points, IQgando et al. (O05) as cyan diamonds. iReda et alj j 2007) (R07) as orange pentagons an d 

|Proctorll2003h (P03) as green triangles. The blue and red arrows are for the results in lWu et al. I d2005f ) and lLa Barbera et al.l (2009). 
Top-right. Our ETGs divided in young and old central ages are compared with literature data: diamonds with continue bars are for ETGs 
with < age\ 6 Gyr; triangle with dashed bars are for agei > 6 Gyr. Bottom-left. Legend of symbols for literature data. Bottom-right. 
Age gradients divided by age are plotted versus literature data (M03, SB07 and R09) as above. 



(as in lKawata fc Gibsonllioolh . but it fails to reproduce our 
observed gradients for the very massive systems. Moreover, 
as seen in Fig. [H the absolute value of Vz of ETGs is lower 
than for the LTGs, probably as a consequence of the dilution 
of the gradient due to the higher-density environment where 
ETGs generally live in (see also discussion in t )3.2[) . It might 
be also the ef fect of other mec hanisms at work such as merg- 
ing (see, e.g., lKobava shi 2004), which however is expected t o 
be less effective in these mass ranges ijde Lucia et al. I 2OO6K 

At very low masses (M* < 10 9 ' 5 Mq and loga c < 
1.8 km/s) the Vz turns to even positive values, which mainly 
corresponds to systems with negative V agc as in Fig. [5] and 
very low central metallicities as in Fig. [6] these dwarf-like 
systems are comp atible with the expanding shell model from 
iMori et all (|l997h (see also Fig.©. 

At larger masses, the shallow metallicity gradients of 



ETGs suggest that these have experienced merging and/or 
tidal interactions (at a rate that could increase as a func- 
tion of the stellar mass) which have diluted the Vz- Such 
events might have taken place in earlier phases of the galaxy 
evolution, as indicated by the presence of null age gradients 
in old systems. In Fig. [S] we compare the V age to the pre- 
diction from the hierarchical simulations in Irlopkins et al.l 
2009a, where a perfect agreement is found for the old sam- 
ple (see also the good agreement with literature data if 
Fig. [7]). From the same simulations we also see that sys- 
tems with younger cores are expected to show positive gra- 
dients which are fully consistent with our low central age 
systems as in Fig. [4] and also reported in Fig. [8] In the 
same figure, metallicity gradients a re in qualitativ e agree - 
ment with gradients of remnants i n iHopkins et all (|2009af) 
and the medians from iKobavashil l|2004l ). while the merg- 
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Figure 8. Left. Age gradients versus simulations. Comparison with gas-rich merging predictions in iHopk ins et al. (2009a). with the 
darker (lighter) region s howing galaxies older (younger) than 6Gyr. Middle. Metallicity gradients versus si mulations. The pred i ctions 
from merging models in iBekki fc Shioval h999ft are shown as short dashed line, dissipa tive collapse models inlKawata fc Gibsonl J2003h 
as long dashed line, the remnants of major-mergings between gas-rich disk galaxies in IHopkins et al.l ]2009aj)_ as yellow shaded region, 
the typical gradient for major (continue gray line) and non-maj or (dashed gray line) mergings in lKobavashil J20 04), The violet thick line 
shows the pred ictions as in the simulation of iMori et al] ljl997T l: blue and pink lines are the result from the chemo-dynamical model in 
iKawatal j200ll ), respectively for strong an d weak supernovae fee dback models (the B-band magnitude in this work is transformed in a 
velocity dispersion cro using the relation in lTortora et al] l|2009al )). The symbols are as in Fig. [Hand the dashed red line is for the whole 
ETGs sample. Right. Symbol legend for simulations. 



ing models in IBekki fc Shioval (|l999l ) reproduce the mean 
gradients in galaxies with log cro J>2.4, but fail at lower ctq. 
Once again, if we consider the older systems only, they are 
in a better agreement with the merging simulations: this 
is expected since, after the initial gas rich-merging events 
that might produce b oth a larger cent r al metallicity and a 
posit ive age gradient |Kobava shi 2004, Mi hos fc Hernquistl 
119941 ). subsequent gas poor- merging may dilute the posi- 
tive age gradient with ti me as well a s make the metallicity 
gradients to flatten ou t i]Whitelll980l , IHopkins et al.ll2009al . 
|Pi Matteo et all2009h - see also left column of Fig. [6] In this 
respect, massive and old ETG systems seem to be fully con- 
sistent with the merging scenario. However, a further mech- 
anism that may act to produce the shallower (or almost null) 
color and metallicity gradients in old massive ETGs, with 
M* ^lO 11 Mq and logcr ^2.4 (Fig. ij), mi ght be due to 
the st rong quasar feedback at high redshift (Tortora et al.l 
l2009bl ). while steeper metallicity gradients at lower masses 
could be linked to less efficient AGNs. 



5 CONCLUSIONS 

We have investigated the colour gradients in a sample of 
~ 50 000 local galaxies from the SDSS as a function of 
structural parameters, luminosity, and stellar mass. CGs 
have been found to correlate mainly with luminosity and 
stellar mass. They have a negative minimum (V g _i ~ 
—0.2) at M* ~ 1O 1O ' 3 M0, and increase with the mass for 
M* ^10 10 ' 3 M Q , the very massive galaxies, Af» ^10 n M Q , 
having the shallower values (~ —0.1). On the other mass 
side, the gradients decrease with mass and become positive 
for M» < 10 9 Mq . These trends are mirrored by similar be- 
haviours with luminosity and galaxy size, e.g. the gradients 
have a minimum at r ~ —20 mag and log R c g ~ 0.5, then 
increase toward the small and the large end of the parameter 
distribution, turning to positive values for log i? c ff <C, and 



r ;> — 18 mag. The dependence on velocity dispersion is very 
loose. A clear dichotomy is also found when looking at the 
dependence on the Sersic index n which suggests a distinct 
behaviour between the ETG and LTG. In fact, these two 
families mark clear differences in their trends with struc- 
tural parameters (e.g. mass and <jq). 

For LTGs, V g _i monotonically decreases with the mass, 
with more massive systems having the lowest CGs (~ —0.4) 
and less massive systems (M* < 1O 8_8 ' 5 M0) showing 
even positive gradients. ETGs have negative gradients 
mildly increasing with mass for logM* ^io 10 - 3 " 10 - 5 Mq, 
which marks the mass scale for the gradient slope in- 
version (see Table [1]). This result is consistent with 
ISpol aor et afl l|2009l ) and reminiscent of the typical mass 



in galaxies drastically change (Capaccioli et al. 1992a: 


Capaccioli. Caon & D'Onofrid Il992bl, 


Kauffmann et al. 1 


20031, Croton et all 200ot Cattaneo et al. 


2008). A similar 



trend is observed when plotting the gradients as a function 
of effective radius, while a tighter trend with velocity 
dispersion is evident for ETGs only (Fig. [4j. This was 
masked when plotting LTGs and ETGs together. 

We have used galaxy colors at 0.1i? e ff and lR e s 
in our synthetic spectral models to determine the vari- 
ation of age and metallicity at these radii. The ob- 
served trends of the CGs with mass and <7o are corre- 
lated with a similar trends for metallicity and age gra- 
dients. Despite the large scatter of the data, the strong 
correlation of metallicity gradients with the central ve- 
locity dispersion of ETGs is clear, with a turnoff point 
at log cro ~ 2.2 km/s (see Table [T]). These results are in 
very goo d agreement with a collection of r e sults from lit- 
eratu r e ijMehlert et al] 120031. iProctorl 120031, lOgando et al] 



20051. lReda etal.1 l200l ISanchez-Blazquez et al. I 120071, 
Koleva et all l2009al lbl ISpolaor et al] 120081 . iRawle et all 
2009), in particular for the old galaxies in our sample. 
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A remarkable result of our analysis is the confirmation 
that the galaxy (central) age is one of the main drivers of the 
scatter of the age and metallicity gradients, with the older 
systems showing generally the shallower gradients with re- 
spect the young ones. In Fig. for instance, we show that 
in the low mass regime, at fixed on, older galaxies have on 
average shallower gr adients, consistently with the results in 
ISpol aor et all (J2009), while systems with late formation or 
with a recent SF episode have a larger spread (and metallic- 
ity gradients do wn to very low value s, ~ —0.6), consistently 
with findings in lKoleva et alj |2009al fbT). The measured scat- 
ter might be the consequence of a variety of phenomena 
affecting the dwarf galaxy evolution, such as (soft) SN feed- 
back, interaction/merging in the high density environment, 
and star- formation by shell expansion (Fig. 8). On the other 
mass side, a further factor of the spread of the massive sys- 
tems gradients might b e the randomness o f the initial con- 
ditions of the mergings (|Rawle et alj|2009h . 

Previous attempts to quantify the correlations of gradi- 
ents with luminosity, mass, or velocity dispersion have often 
failed, mainly because of the exiguity of the galaxy sam- 
ple (e.g. iPeletier et al.l Il990al . iKobavashi fc Arimotol 1 19991 . 

fc Ohtall2003| )l Only recently a clear correlation of 



Ta 



the color and metal licity gradients wit h mass has been as- 
certained (see e.g., iForbes et al.l [2005) . poin ting to differ- 



ent t r ends for high and low mass galaxies l|Spolaor et al.l 
2009, iRawle et alj|2009h . Our analysis has confirmed and 
reinforced these results, as it relies on one of the largest 
local galaxy samples including dwarf, normal, and giant 
galaxies. Using CGs derived by the structural parame- 
ters in B05, we obtained statistically meaningful trends, 
which allow us to fix the link with physical phenomena 
as well as to ma ke direct compar i sons with predictions 



from simulations ilMori et al.l Il997l iBekki 
Kawatal l200ll iKawata fc Gibson! 120031 
Hopkins et afll2009al ). ~ 



Shiova 

Kobavashil 



1999 



2004 



We can finally draw the physical scenario sketched in 
Fig. [9] resting on the results discussed so far. The forma- 
tion of LTGs and less massive ETGs is mainly driven by a 
monolithic-like collapse as they have almost null or positive 
age gradients and negative metallicity gradients decreasing 
with the mass. These are in fact well explained by simple 
gas inflow and feedback from SN and evolved stars, and 
also reproduced in simulation s of dissipative col l apse and 
SN feedback models in Fig. [8l (lLarsonlll974 Il975l, ICarlberd 
11984 lArimoto fc Yoshiilll987UKawata fc Gibsonll2003l ). The 
difference in the magnitude of the gradients between LTGs 
and lower mass ETGs mainly resides in the effect of the envi- 
ronment as the former live in very low density environments, 
while the latter stay in massive haloes, together with more 
massive ETGs. Here they experience strong gravitational 
interactions producing shallower color and metallicity gra- 
dients. The efficiency of such phenomena is strong enough 
for lower mass ETGs, th us shaping the observed d ecreasing 
trend with mass (see e.g. IDekel fc Birnboimll200^ ). 

The most massive ETG systems have shallower gradi- 
ents ilBekki fc Shiovalll999l . lKobavashill2004 iHopkins etall 
12009a ). flattening out with the mass due to the in- 
creasing intervention of (gas-rich and -poor) merg- 
ing, tidal interactions, an d quas ar /radio mode AGN 
feedb ac k dCapaccioli et al.1 [l992al. iFerrarese fc Merrittl 
|2000| . Ide Lucia et al. I I2006L IDekel fc Birnboiml 120061 . 
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Figure 9. Color-mass diagram. The contours show the density 
of data-points with gray scale going from darker (low density of 
galaxies) to brighter regions (high density). The thin white line 
sets the separation between RS and BC; the thick vertical one 
gives the mass scale, M, ~ 10 10,3 , which separates galaxies be- 
longing to the RS in normal and dwarf ETGs, and sets a quali- 
tative upper mass for LTGs. The arrows give information about 
the efficiency of the phenomena which drive the two-fold trend 
we discuss in the paper. AGN feedback and merging are impor- 
tant at high mass with an efficiency increasing with mass, while 
SN feedback and gas inflow drive the galaxy evolution in the less 
massive side, with an efficiency that is larger at the lowest masses. 
Galaxies in RS and BC lie in environments with a different den- 
sity, which is manifested as a difference of the gradients in the 
two samples. 



Liu. Jiang, fc Gul 120061 . ISiiacki et ail 120071 . ICattaneo et~al] 
2008? ) . All the simulations collected in literature fail to 
reproduce the fine details of the trends we find, suggesting 
that a full understanding of physical processes involved in 
the galaxy evolution is still missing. 

In future analysis we plan to enlarge the wavelength 
baseline and also to use line-strength measurements for the 
stellar model (when available). We will investigate the sys- 
tematics induced when other synthetic prescriptions are as- 
sumed and the dependence of the CGs on the galaxy star 
formation history. We will derive the M/L gradients, if any, 
and discuss their impact on the dark matter content of ETGs 
and correlate the CGs with the DM fraction of these systems 
as probe of the galaxy potential wells. Finally, the realiza- 
tion of hydrodynamical simulations of jets emitted by AGN 
would be useful to shape, together with galaxy merging, the 
color and metallicity gradients for massive galaxies. 
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APPENDIX A: SAMPLE SYSTEMATICS AND 
CONTAMINANTS 

We analyze various systematics in the galaxy sample and in 
the selection criteria. 



Al Sample systematics 

We quantify the effect of source incompleteness on our re- 
sults, in particular on the age and metallicity gradients as a 
function of stellar mass. 

Our sample is magnitude limited, and misses increas- 
ingly fainter and less massive galaxies at larger redshifts; 
at z sr>ec = 0.02, the lowest mass galaxies have M* ~ 



e.g. di l ^spec 

10 8,2 Mo, while at z avec = 0.05 this limit increases 
M* ~ 10 9 M Q (see lBlanton et al.|[2005al and the left panel 
in Fig. IA1|) . We already mentioned that dwarf galaxies can 
be missed due to surface brightness selection effects (pho- 
tometric incompleteness). Other less obvious sources of in- 
completeness are defects in deblending and the absence of 
reliable redshift measurements. 

By analyzing the s urface brightn e ss com pleteness of 
this sample of galaxies, Bl anton et al.l (|2005af ) found that 
for r < — 18 mag (M* < 1O 9 M ) the catalog is com- 
plete at more than 90%, while at fainter magnitudes the 
photometric pipeline could mistakenly deblend galaxies or 
overestimate the background sky level. However, down to 
M« ~ 10 8 5 Mp the overall incompleteness is not dramatic 
jBlanton et al.l l2005al . iBaldrv. Glazebrook fc Driveij l2008t 
iLi fc White! 120091 ). Moreover, at z apec < 0.02 the distribu- 
tion of galaxies seems patchy and the fraction of massive 
galaxies with M* :>10 9 Mq turns out to be smaller, and not 
homogeneously distributed as the one at larger redshifts. In 
order to account for all possible biases due to these miss- 
ing galaxies we have redone the age and metallicity gradient 
plots for the subsamples of galaxies having z apec > 0.02 
and Zspec > 0.03, in order to avoid the regions with higher 
incompleteness. The trends found for the two subsamples 
are almost identical to the ones for the full data set, which 
suggests that the impact on our results of all sources of in- 
completeness is negligible, both at low and high masses. 
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Figure A2. Systematics in the V age and Vz as a function of 
stellar mass due to changes in selection criteria. Boxes and bars 
are for our reference results. Continue, long-dashed, and short 
dashed lines are for a) C > 2.6, b) C > 3 and c) 2.5 < n < 5.5 
respectively, and triangles for our reference criteria with a cut in 
the total color (g — r > 0.55). 



It has been suggested that the Sersic fitting procedure 
adopted in B05 contain significant biases. In particular, at 
very high luminosities, Sersic indices, effective radii and 
fluxes can be underestimated by ~ 0.5, ~ 10 — 15% and 
~ 10%, respectively (B05). Similar considerations hold for 
systems with high effective radii and Sersic indices. The 
Sersic p rocedure u s ed in B05 have been extensively ana- 
lyzed in iGuo et al.l l|2009l ). who have shown that the main 
source of systematics is introduced by the ID Sersic fit and 
background sky level estimate. However, these systematics 
turn out to produce shallower CGs (of ~ 0.03) than our es- 
timates, thus leaving the main trends with luminosity and 
mass unaffected. On the other hand, for the systems with 
very low Sersic indices, the B05 fitting procedure recovers 
quite well the parameters, although a slight under estimate of 
fluxe s and an overestimate of n might still occur (|Guo et alj 
120091 ). 



A2 Morphological separation and contamination 

Seeking for a trustful morphological classification, the eye- 
ball check or the B/T are the best approaches, but un- 
fortunately they ar e not suitab l e procedure for very larg e 
sample of galaxies (|Allen et all |2006| . iLintott et al"1l2008h . 
Indeed, the use of some objective automated morphologi- 
cal classification is more efficient. For these reasons, we de- 
cided to adopt a classification with the minimal parame- 
ter set to minimize the probability of erroneously discard- 
ing ETGs, even if this would make the LTG contamina- 
tion non-negligible. We have "a posteriori" checked the ef- 
fect of such a contamination on the results. We have ver- 
ified that the use of the concentration and Sersic n— index 
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Figure Al. Left Panel. Spectroscopic redshift as a function of stellar mass for the B05 galaxy sample. The two lines are relative to 
the redshift thresholds we use to select subsamples of galaxies. Middle Panel. Age gradients as a function of stellar mass. Right Panel. 
Metallicity gradients as a function of stellar mass. The red and blue symbols with bars are relative to the full galaxy sample for ETGs 
and LTGs respectively, as discussed in the text, and the continue and dashed lines are for subsamples with z sp ec > 0.02 and z sp ec > 0.03. 



parameters is efficient to this purpose. In fact, the Sersic 
index is a measure of the steepness of the surface bright- 
ness profile and is historically considered as a good indi- 
cator of galaxy type: LTGs are well fitted by light profiles 
with lower n, being n = 1 the prototype of light profile 
for galaxy disks, and n > 1 the ones typical of ETGs. Re- 
cently, the concentration parameter determined from the 
SDSS pipeline has be en claimed to be an even better proxy of 
galaxy morphology ( Shimas aku et al.l l2~00ll. IStrateva et al.l 



1200 1| . iNakamura et al. 1 120031 . IWeinmann et al. 1120091) . We 
have checked that the results on the gradient trends do 
not strongly depend on the constraints adopted on the 
two parameters for the ETG/LTG separation (see Fig. IA2|) 
with some significant change raising for age gradients at 
Af» < 10 9 ' 5 Mq. We also analyzed the impact on our re- 
sults when a cut in colour is applied. 

As anticipated, our selection criteria are not meant to be 
contamination free, with ETGs being possibly contaminated 
by disk galaxies: e.g., both n and C depend on the inclina- 
tion and edge-on disk or spiral galaxies can show parameter 
values typical of ETGs. An approach to avoid this kind of 
contamination can rely on the axial ratio, but i n this case we 
woul d loose edge-on SOs misc lassified as LTGs (jMaller et al.l 
2009, Wei nmann et al. 1 ^009) . In order to evaluate the frac- 
tion of contaminants in our ETG sample we have made an 
eyeball inspection of a subsample of 300 randomly selected 
galaxies. This classification may be still prone to subjectiv- 
ity, thus four of us have checked the galaxies and produced 
an average classification. Eight galaxy images are shown in 
Fig. El We find that ~ 70 - 75% of galaxies are ETGs, 
while the remaining galaxies are uncertain objects or clear 
spirals. These uncertain or misclassified galaxies mainly fill 
the intermediate-mass region and many of them have red col- 
ors g — i J>0.5. Although the fraction of such contaminants 
is not negligible, we have found that the gradients are not 
strongly biased on average, as shown in Fig. IA4I The bot- 
tom line of this check is that the inclusion of contaminants 
leaves the average gradients unaffected. For instance, on the 
mass range 10 < logM, < 11, we obtain as median value 
Vago = 0.12, and 0.03 for the ETGs, and the contaminants 
respectively, as shown in Fig. IA4I while we find Vz = —0.47 
and —0.54 for the same systems. However, when comput- 
ing the median values for the whole sample, the result is 
unchanged from the one of systems correctly classified as 
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Figure A5. Age (left panel) and metallicity (right panel) gradi- 
ents as a function of stellar mass for ETGs (red symbols), LTGs 
(blue symbols) and ITGs (green symbols). 



ETGs. No net difference in the central galaxy age is found, 
while central metallicity for contaminants is larger than the 
one for ETGs, consistently with the steeper metallicity gra- 
dients for galaxies with larger metallicity (see Fig. [6} . We 
have also checked that the presence of contaminants in ETG 
sample, in particular, does not affect the dependence on the 
central age of V ag e and Vz as discussed in Fig. [4] Thus, we 
can safely assume that the central age is a genuine param- 
eter which rules the scatter of the relations with mass and 

(TO- 

Finally, the conservative choice of both n and C pro- 
vided a quite sharp separation between ETGs and LTGs. 
For sake of completeness, the residual galaxy sample with 
a mixed morphology, i.e. ITGs, turned out to have also an 
intermediate distribution of the age gradients as shown in 
Fig. IA5I On the contrary, metallicity gradients of ITGs have 
a trend which is more similar to the one of ETGs. 
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Figure A3. Images of 8 galaxies classified by eyeball checking. The galaxies in the top panels are classified as E/SO, the first one in the 
bottom panels is uncertain, while the last three are clear example of spirals. 
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Figure A4. Results for galaxies extracted for eyeball check. The red points are for galaxies classified as ellipticals, while blue ones are 
for the galaxies with an uncertain classification or spirals. Left Panel. Color-mass diagram. Middle Panel. V age as a function of stellar 
mass. Right Panel. Vz as a function of stellar mass. 



APPENDIX B: SYSTEM ATICS IN STELLAR 
FITTING 

We examine the role of systematic uncertainties on the stel- 
lar populations results. First, we address the problem of de- 
generacy in synthetic spectral models (§ IB1|) . then derive 
the optical-infrared and infrared CGs from our fitted stel- 
lar populations and check the consistency with results from 
literature (S IB2I) . Dust extinction is discussed in S IB3I 



Bl Parameter degeneracies and wavelength 
coverage 

We have analyzed the optical CGs, fitting the internal and 
external colors with synthetic spectral models to derive 
the age and metallicity of stellar populati ons. Because of 
the well known age-meta l licity degenerac y ( Wort hey) Il994l . 
iBruzual fc Chariot IlioOol iGallazzi et al. 1120051 ). the stellar 



parameters and the relative gradients might be severely bi- 
ased, in particular by using the optical colors only. Increas- 
ing the wavelength coverage or adding spectral information 
would be the best approach to the problem. 

Near-IR bands are primarily sensitive to the red stellar 
populations which represent the main fraction of stars form- 
ing an evolved galaxy, allowing to trace the visible mass of 
old galaxies, and are less affected by dust extinction. Their 
inclusi on is suitable t o alleviate the age-metallicity degen- 
eracy (Ide Jond 1996J, ICardiel et all 12001 iMacArthur et al.l 
|2004| . [Chang et ail [2006'). but is not expected to be defini- 
tive. Different colors a re found to hav e a different sensitivity 
to age and metallicity (|Li et al.ll200^ . For instance, both op- 
tical and infrared colors are found to be sensitiv e to metallic- 
ity, while age mostly a ffects the optical colors (|Chang et alj 
l200rj) . |Wu et all (|2005l ) showed that the inclusion of near-IR 
bands allows one to recover much more accurate stellar pa- 
rameters from the fit, but the confidence contours are still 
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Table Bl. Results from Montccarlo simulations. We show the 
median relative age A(age) = (agefu — agei n )/agei n and metal- 
licity A(Z) = (Zfn — Zi n )/Zi„, where fit and in are for estimated 
and input parameters. The results for different initial perturba- 
tions 5 and without or with IR photometry are shown. 
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Figure Bl. Comparison of age (left) and metallicity (right) gra- 
dients estimated from Montecarlo repetitions (y-axis), given the 
initial input values (x-axis). Here we consider our extreme error 
budget, 8 = 0.05. Red and dashed blue contours are the isoden- 
sity corresponding to the results from the stellar population fit, 
using only optical and optical+infrared colors respectively. 



affected by the age-metallicity relation, and the gain with 
respect the optical bands is minimal. It is reasonable to ar- 
gue that the use of a large statistical sample, as we adopt 
in this paper, is more effective in reducing the overall un- 
certainties of the galaxy properties than the adoption of a 
larger filter baseline. 

Of course, the adoption of spectroscopic line absorp- 
tion indices would be the best indicators of stellar popu- 
lation parameters and gradients, as they permit to break 



the age-metallicity degeneracy (e.g., [B urstcin c t al. 19841. 
Henry fc Worthevl 1 19991. ICarollo et all 1 19931 iDavies et all 



199.i l Mehlert et all B . However. sne c^Tin dices of large 
galaxy sample are still difficult to obtain for the observa- 
tional limits imposed by the low starlight fluxes (e.g. around 
R e fi and beyond). Recently, integral-field spectroscopy has 
been used to determine gradients in line absor ption indices 
(e.g. iKuntschner et aT I l2006l . iRawle et al. I l2009h of elliptical 
and lenticular galaxies, but this approach is only applicable 
to limited samples of systems. 

Despite the limited wavelength coverage, but taking ad- 
vantage of the large statistical sample we have shown that 
both spectroscopy and photometric analysis converge to sim- 
ilar result s for stellar popula tion gradients (see e.g. £|3. 4[) . 
Similarly, Rawle et al. (2009) have shown that optical col- 
ors and CGs derived from spectroscopic age and metallicity 
are generally consistent with values directly recovered from 
photometry. Moreover, the recovered correlations between 
metallicity and age gradients are not spurious and not gen- 
erated by correlated measurement errors on these quantities. 

However, in the attempt of putting our result further 



on a more solid ground, we have tested the reliability of 
our color modelling technique and checked for the presence 
of spuriously generated correlations between age and Z or 
gradients V agc and Vz , running several Monte Carlo simula- 
tions. We extracted 1000 simulated galaxy spectra from our 
BC03 SED libraries with random (agei, Zi) for i — 1,2 (i.e. 
with no correlations among these parameters and with Vz , 
and Vagc randomly distributed) , and applied our fitting pro- 
cedure, comparing the estimated parameters with the input 
model values. We perform the fit 1) only using the optical 
SDSS bands ugr iz and 2) adding th e JHKs photometry from 
2MASS survey l|jarrett et al.|[2003T ). to understand the sys- 
tematics that the fit of optical colors can induce. Typical 
observational uncertainties for each band are assigned (i.e. 
Su ~ 0.06, Sg ~ 0.035, Sr ~ 0.035, Si ~ 0.035, Sg ~ 0.04, 
8 J ~ 0.08, 8H ~ 0.12, 8Ks ~ 0.09, from matched catalog in 
B05) and the input model colors are perturbed adding a ran- 
dom step in the interval (—5, +8) with 8 — 0.01, 0.03, 0.05 to 
account for measurement systematics. Results of this analy- 
sis are collected in Table iBll and Fig. IB1I We found that, on 
average, both the stellar parameters and gradients are fairly 
well recovered from the fitting procedure, with a scatter in- 
creasing with 8, and no spurious correlatio ns are induced at 
mor e than 99% confidenc e level (see also IWu et al"] (|2005l ) 
and lTortora et alj (l2009ah for similar analysis). We confirm 
the finding of IWu et al. 

1 120051 ) 

that the inclusion of near-IR 
bands allows a better recovery of the input stellar parame- 
ters with a scatter in the range 5 — 25%, while the scatter 
is larger (in the range 10 — 50%) when only SDSS bands 
are used. Thus, the most remarkable result of this analysis, 
is that, notwithstanding the larger scatter, we are able to 
successfully recover the stellar population parameters, even 
without near-IR photometry. 

B2 Comparison with IR CGs 

One of the advantages of synthetic spectral modelling is to 
allow the recover of synthetic IR colors by using the best 
fitted spectral parameters. Therefore, we derive the colors 
g — J , g~H and g — Ks at 0.1i? c ff and 1R c b and the relative 
CGs, which show trends with mass similar to those discussed 
in this paper. If we select luminous (r < —20 mag) ETGs 



-0.33, V g - H 



-0.37 



and V„-Ks 



-0.40, with 25-75 percentiles in sample dis- 



tribution of 0.17, 0.18, and 0.20, respectively. As already 
extensively discussed in the text, to com pare these syn- 
thetic CGs with the observed values in lLa Barbera et akl 
(2009), we select older systems with agei > 6Gyr. In 



this case we have V 



g-J 



-0.26, Vg-H 



-0.29, and 



Vg-Ks = -0.31 (with scatter of 0.16, 0.18, and 0.20), which 
are fairly consistent with the results reported by these au- 
thors. A fairly good agree ment is also found with the median 
Vg-Ks = -0.29 ± 0.07 in lWu et all (|2005T l. 



B3 Internal dust extinction 

In our analysis we have assumed no internal extinction. 
This assumption is more reasonable for ETGs, where only 
a negligible f raction of interstellar matter is in the form 
of dust (e.g. Kormcndy fc Diorgovskil [19891 , IWise fc Silval 
1 19961 . [5 orman et al l l2003l l7 but it might be unappropri- 
ate for LTGs, as these systems have a larger amount 
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Figure B2. Gradients for the spectral model leaving free the galaxy age, metallicity, and dust content. The open symbols with bars are 
for our reference fit, while the line is for the fit with dust extinction free to vary. Left Panel. V ag e as a function of stellar mass. Middle 
Panel. Vz as a function of stellar mass. Right Panel. Gradient in dust extinction as a function of stellar mass. 



of internal dust (e.g. iBrammer et al . 2009). Usually no 
dust extinction is included in the analysis of ETGs 
llTamura et al-lfeoOOllTaniura fc Ohtall2000l . rTamura fe Ohtal 



I200&1 Wu et alj l2005h. but dust gradients wou ld have 



effects in CGs (e.g. iGoudfrooii fc de Jond Il995h . How- 
ever, for ETGs it has been shown that metal absorp- 
tion lines gradients, b eing less affected by dust extinc- 
tion than colors (e.g. iMacArthurl l2005t ). are unlikely to 
be produced by pure dust extinction gradients to mimic 
the population gradients inferred by broadband colors 
jPeletier et aljl990al.lpavies et al1ll993l.lCarollo et aljl993t 
iKobavashi fc Arimotolll999l iRawle et alj|2009h . Thus they 
are not dominant but could still contribute to a fraction of 
CGs. Similarly, in LTGs the dust seems to have a n egligi- 
ble effect in shaping CGs (e.g. iMacArthur et al1l2004l). even 
thoug ht there are still open questions ( MacArthur et all 
2009) which might lead to a different conclusion. To ana- 
lyze the impact of dust extinction on our results, we fit the 
colors at 0.1R c b and li? c ft using a library of spectra with 
age, metallicity and dust content free to change. The ex- 
tinction curve of ICardelli. Clayton fc Mathisl i| 19891 ) is used 
to account for the extinction contribution in different wave- 
bands and the color excess E(B — V) in the range — 0.5 is 
used as a free parameter. The results are shown in Fig. IB2I 
where they are compared with our reference results with no 
dust. The main trends are unaffected by the dust, although 
we note some differences in the V ag o at M* J>10 9 ' 5 Mq. 
On the other hand, only slight differences are observed for 
metallicity gradients, and the qualitative trends are unaf- 
fected both in ETGs and LTGs. Thus this analysis has 
shown that, although dust extinction is not negligible in 
both ETGs and LTGs, the trends and the amplitude of gra- 
dients are still dominated by changes in stellar population 
parameters. We show the dust extinction gradients, defined 
as Ve(b-v) = E(B - V) 2 - E(B - V)i in the right panel 
of Fig. IB2I On average, the ETGs exhibit null or slightly 
positive dust gradients (with a tail towards negative values, 
too) , while LTGs have opposite behaviour with a larger dust 
content in the inner regions, showin g steeper dust gradients 
at intermediate high masses (e.g., iDalcanton fc Bernstein! 
2002. IMacArthur et alJ fcOQ^. Finally, we have checked the 
absence of any correlation between age and metallicity gradi- 
ent with both V e(b-v) an d the central extinction for LTGs 
which has suggested that the results of these systems are not 



affected by the presence of dust. We remark that the use of 
only optical bands and the addition of IR photometry is not 
sufficient to break the degeneracy between stellar popula- 
tion parameters and dust, for this reason we are tempted to 
consider these results more as a caveat for this particular 
data-set rather than a conclusive check. If the correction of 
the age gradients for the dust extinction is right, the net ef- 
fect is that the age gradient might be less affecting the color 
gradient for massive systems, and only metallicity gradient 
should matter. 



